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Abstract. The description of electron-electron interactions in transport problems 
is both analytically and numerically difficult. Here we show that a much simpler 
description of electron transport in the presence of interactions can be achieved 
in nanoscale systems. In particular, we show that the electron flow in nanoscale 
conductors can be described by Navier-Stokes type of equations with an effective 
electron viscosity, i.e., on a par with the dynamics of a viscous and compressible 
classical fluid. By using this hydrodynamic approach we derive the conditions for the 
transition from laminar to turbulent flow in nanoscale systems and discuss possible 
experimental tests of our predictions. 



1. Introduction 

The electron liquid is both viscous and compressible; properties which suggest an 
intriguing analogy with a classical liquid [lj. This analogy is even more compelling 
when one recalls that the time-dependent many-body Schrodinger equation (TDSE) 
can be cast, quite generally, in a "hydrodynamic" form in terms of the density and 
expectation value of the velocity operator [21 E]. In the classical case one can derive 
time-dependent equations, called Navier-Stokes equations, for the velocity field of the 
fluid as a function of its density, visco-elastic coefficients, pressure and the geometric 
confinement [lj. These equations are centrefold in hydrodynamics and describe both 
laminar and turbulent regimes. If we could derive similar equations in the quantum 
case we would have a powerful tool to investigate a plethora of effects related to 
electron-electron interactions on a much simpler level than solving for the many-body 
Schrodinger equation. Unfortunately, the derivation of these equations in the quantum 
case is generally not possible. 

In this paper we show that transport in nanoscale conductors satisfies the conditions 
to derive quantum Navier-Stokes equations. This is simply due to the geometric 
constriction experienced by electrons flowing in a nanostructure which gives rise to 
very fast "collisional" processes [HE]. In this regime, we show that one can truncate 
the infinite hierarchy of equations of motion for the electron stress tensor to second 
order and thus derive quantum hydrodynamic equations. With these equations we 
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first rederive conductance quantization in a quasi-lD non-viscous, incompressible fluid 
thus making the connection between quantum transport and this hydrodynamic picture 
clearer. We then predict the conditions for the transition from laminar to turbulent 
flow in quantum point contacts (QPCs) and suggest specific experiments to verify our 
predictions. 

Let us start from a general many-body Hamiltonian H = T + W + V ex t, where T is 
the kinetic term, V ext an external potential and 

W = — J dr J dr'ip^ \r)ip\r')w{\r — r'\)ip(r')ip(r), (1) 

where ip(r) are field operators and w(\r — r'\) is the Coulomb interaction potential. It 
is well known that the TDSE can be equivalently written as two coupled equations of 
motion for the single-particle density, n(r, t), and velocity field, v(r, t), as obtained from 
the Heisenberg equation of motion for the corresponding operators [2j [3]. The single 
particle density is defined in terms of the field operator as n(r,t) = (i[)'(r,t)i/;(r,t)), 
while the velocity is given by v(r,t) = j(r,t)/n(r,t) where the current density is given 
by j(r,t) = eh{[il)'(r,t)ViJ)(r,t) — V"0'(r, t)ip(r, t)])/2mi. For clarity we rewrite the 
Heisenberg equations of motion here (summation over repeated indexes is understood) 

D t n(r, t) + n(r, t)Vv(r, t) = (2) 
mn(r, t)D tVj (r, t) = -V.P^r, t) - n(r, t)VjV ext (r, t) (3) 

where D t = dt + v ■ V is the convective derivative, m the electron mass, and Pij 
= rnj) = x,y,z) is a stress tensor, exactly given by the sum of kinetic and 
interaction stress tensors [21 E] - The interaction tensor is 

WiAr,t) = -~ [ dr'^^jP- f 1 dX G 2 (r + Xr', r - (1 - A)r') (4) 
2 j | r /| «9| r /| j Q 

where G2(r,r') = (%p\r)h(r')ip(r)) is the two-particle density matrix and A a parameter 
that defines the geodesic which connects two interacting particles [3]. As mentioned 
above, ([3]) has an appealing "hydrodynamic" form where all many-body information 
is included in the stress tensor Pjj. In order to solve ([3]) one proceeds by calculating 
an equation of motion for which can be derived from the Heisenberg equation of 
motion of the particle creation and destruction operators. However, this equation of 
motion contains the three-particle density matrix. In turn, the equation of motion for 
the three-particle density matrix contains the four-particle density matrix and so forth, 
thus generating an infinite hierarchy of nested equations, making the problem practically 
unsolvable [21 [3] . 

2. Quantum Navier-Stokes equations 

We show here that in the case of electrical transport in nanoscale systems we can instead 
close this set of equations. We proceed as follows. First, let us derive the dependence of 
the stress tensor P it j on the rate at which the system reaches a quasi-steady state. We 
can obtain this dependence from the quantum kinetic equations for the nonequilibrium 



Hydrodynamic approach to transport and turbulence in nanoscale conductors 



3 



distribution function f(r,p,t) (p is the momentum), which can be derived from the 
TDSE equation with standard techniques, [fjjij] in a co-moving (Lagrangian) reference 
frame moving with velocity v(r,t) [7] 

/[/] = D t f(r, P , t) + >. Vf(r, P , () + eV^^M - p • V„^M 
m op op 

_ mD /ttMA (5) 
op 

where / is the usual collisional integral [6j, ip is the sum of the external potential and the 
Hartree part of the interaction potential. The collisional integral contains two terms, 
one elastic and the other inelastic. In what follows, it is important to realize that both 
terms can drive the system toward a local equilibrium configuration. 

The first two moments of the distribution give the density n(r,t) = J2 P f(.Vi r it) 
and the condition ^2 p pf(p,r,t) = 0. The two-particle stress tensor is related to the 
distribution function as Pij = ^2 p PiPjf{p,r,t)/m. Similarly, higher order moments of 
the distribution function produce higher order stress tensors. Introducing these moments 
into ([5]) and comparing with (J3j) we can write the stress tensor in terms of the collisional 
integral I according to the equation of motion 

— dp I[f}PiPj = D t P iyj + P M V ■ v + P^kVkVj + P^Vm + VkP-%, (6) 
ui j 

where P^ is the three-particle stress tensor. By writing the equation of motion for P^ 
we would again get an infinite hierarchy of equations. It is interesting to point out that 
the theorems of time-dependent density-functional theory establish that the stress tensor 
Pi j is a universal functional of the velocity and density only (see for example, [8] or [9]). 
This implies that the hierarchy of equations for the moments of the distribution function 
can be formally closed to all orders in the electron-electron interaction. However, we 
note that enters in §6§ only through its spatial derivative. If the latter is small 
then the hierarchy can be truncated [TJ. From ([6]) we easily see that this derivative 
is small compared to the other terms whenever 7 = u/(Lmax(w,i/ c )) <C 1. Here u 
is the average electron velocity, L is the length of inhomogeneities of the liquid that 
give rise to scattering among three particles, u the system proper frequency and u c the 
collision rate. The parameter 1/L enters through the spatial derivative of P^ 3 \ to from 
the frequency dependence of the interactions (in the DC limit of interest here u — ► 0), 
v c through the collisional integral /[/] oc —v c {f — fo), where /o is the equilibrium Fermi 
distribution. This derivative is indeed small for transport in nanostructures: When 
electrons move into a nano junction they adapt to the given junction geometry at a fast 
rate, and produce a quasi-steady state and local equilibrium distributions even in the 
absence of electron interactions [U [5]. This "relaxation" mechanism occurs roughly at 
a rate v c = (At) -1 ~ (h/AE)^ 1 , where AE is the typical energy spacing of lateral 
modes in the junction. For a nanojunction of width w we have AE ~ 7r 2 h 2 /mw 2 and 

I Clearly, for the definition of local equilibrium distribution to be valid any length scale entering the 
problem has to be larger than the system Fermi wavelength. 
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At ~ mw 2 /ir 2 h. If w = 1 nm, v c is of the order of 10 15 Hz, i.e., orders of magnitude 
faster than typical electron-electron or electron-phonon scattering rates. The condition 
7 = uj (Lmax(w, u c )) <C 1 thus requires the length of inhomogeneities L ^> 1 nm, which 
is easily satisfied in nanostructures. Note instead that in mesoscopic structures this 
condition is not necessarily satisfied. In that case, the dominant relaxation rate v c is 
given by inelastic effects, i.e. it is of the order of THz, so that for typical lengths of 
mesoscopic systems, 7 ~ 1 in the DC limit. Nonetheless, the above condition could still 
be valid for high-frequency excitations, like plasmons, and/or very low densities, so that 
stress tensors of order higher than two are negligible. 
Neglecting V k P>2k in © 

we can thus derive a form for P^j. Let us write quite 
generally the stress tensor P it j as P i; j = 6ijP — tt^, where the diagonal part gives the 
pressure of the liquid, and tt^ is a traceless tensor that describes the shear effect on the 
liquid. From ([6]) we thus find that the tensor tt^j can be written as (in d dimensions, 



where 77 is a real coefficient that is a functional of the density [7\. We point out that 
([7]) is in fact a particular case of a general stress tensor with memory effects taken into 
account [TOl [TT1 fT2] . In our derivation this is the first non-trivial term of an expansion 
of the stress tensor in terms of the density and velocity field. Consequently the Navier- 
Stokes stress tensor in ([7]) can be seen as the first-order (non-trivial) contribution to the 
exact stress tensor of the electron liquid (see also [3j [12], [TTj). 

Using this stress tensor we finally get from ([3]) the generalized Navier-Stokes 
equations for the electron liquid in nanoscale systems 

D t n(r, t) = — n(r, t)V • v(r, t), 

mn(r, t)D tVi (r, t) = -V;P(r, t) + V i7 r M (r, t) - n(r, t)V t V ext {r, t). (8) 

Equations ([8]) are formally equivalent to their classical counterpart [1] and thus describe 
also nonlinear solutions, i.e. the possibility to obtain turbulence of the electron liquid 
in its normal state. In the examples that follow, we will consider only the case in which 
the liquid is incompressible so that the viscoelastic coefficients are spatially uniform: 
This approximation is practically satisfied in metallic QPCs but needs to be relaxed in 
the case of QPCs with organic/metallic interfaces (see e.g. |13j). In addition, for this 
case the Hartree potential is constant and its spatial derivative is thus zero. Therefore, 
flU) reduce to the Navier-Stokes equations for the density and velocity of a viscous but 
incompressible electron liquid 



d> 1) 




(7) 



D t n(r,t) = 0, 
V • v(r,t) = 0, 
mn(r, t)D t Vi(r, t) 



ViP(r, t) + r/V 2 ^(r, t) - n{r, t)ViV ext {r, t). 



(9) 
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3. Conductance quantization from hydrodynamics 

Let us first show that we can derive from ([9]) the quantized conductance of an ideal 
(r) = 0) quasi-lD liquid. We consider the electron liquid adiabatically connected to two 
reservoirs, and we call vl(r) and (jll(r) the velocity and chemical potential, respectively, 
in the left (right) reservoir, with hl~Hr = eVuas- From (jHJ) we then derive the Bernoulli's 
equation that states the conservation of energy 

| + Al + ^ = !f + + ^ (10) 

2 ml m 



where h^R) is the enthalpy of the left (right) lead4§| . Since we assume the fluid is 
incompressible hi = ho. By defining the flow velocity v = (vr + Vl)/2 and the co- 
moving Fermi velocity!]]] Vp = (vl — vr)/2 we get from ( flOj) the relation 2mvvF = eVuas- 
By definition, the current is given by J = env so that, by using the ID density of 
states, I = emvvp/irh = e 2 Vbi as /h, which, in the linear regime, gives the quantized 
conductance (per spin) Go = I jVuas = e 2 /h. If we assume that only a fraction T of 
electrons is transmitted due to the presence of a barrier in the liquid, we can argue 
that, in linear response, the current is an equal fraction of the current in the absence 
of the barrier, i.e. / = envT. The conductance is thus G = Te 2 /h in accordance 
with the Landauer two-terminal result [H]. Corrections to this conductance in the 
presence of viscosity in 3D have been estimated in [13] and were found to depend non- 
linearly on the gradient of the electron density. Corrections to this conductance have 
been experimentally observed and justified with a hydrodynamical model where the 
conductance depends on the physical properties of the electron flow [T5] . 



4. Turbulence 



We know that the time-independent solutions of the Navier-Stokes equations (jHJ) or (J9j) 
can describe many different regimes, with the non-linear (turbulent) regime generally 
favoured with respect to the laminar one. In fluid mechanics, in order to identify these 
regimes, it is customary to define a key quantity, the Reynolds number Re, as the only 
non-dimensional quantity that can be constructed out of the physical parameters of the 
system, like the density, the viscosity rj, etc. pp. In the quantum case we follow a similar 
convention. For instance in 2D we define Re as Re = Q/r] = ml/erj where Q is the 
total mass current, and / is the average total electrical current. For small Re the stable 
(and stationary) flow is usually laminar while for large Re the flow is turbulent [TJ. In 
the latter case, one should then observe a local velocity field which varies in space in an 
irregular way, and whose pattern is very sensitive to initial conditions. 

§ To derive (fT0|) one makes use of the relation wV • v = Vv 2 /2 — v x (V x v) and by projecting the 
equation of motion on the tangent to the current flow. The enthalpy is defined as h = P/ n pQ . In the 
ID case VxdeO and h(n) — (irhn) 2 /2m 2 . 

|| vp is the Fermi velocity in the reference frame moving with velocity (vr + vp)/2. Obviously the 
Bernoulli's equation is invariant under any Galilean transformation. 
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Adiabatic QPCs - Let us apply these concepts to the transition between laminar 
and turbulent flow in QPCs. The microscopic geometry of these structures is quite 
complicated so that analytical solutions to (IHl) cannot be generally found. However, 
for an adiabatic 2D constriction, the self-consistent confining electron potential which 
enters can be approximated with the simple form y = ky/\x\ — 8 where k > is a 
parameter that controls the rate at which the constriction opens up and 5 > is the 
opening of the constriction [see inset (a) of Fig. [1] . We assume electrons either originate 
from, or enter into, the region {y = 0,x = [0,5]} of that potential to resurface on the 
opposite side of the structure. We know from classical hydrodynamics that, for any 
constriction, turbulent flow may exist only on the side where the constriction acts as a 
source of electrons pQ. We will therefore discuss only this side. We are interested here 
in the case of k large, i.e. true adiabatic limit. Since we are dealing with a viscous 
liquid we assume a no-slip condition at the boundary, i.e., at the boundary the velocity 
component parallel to the boundary is zero. This condition can be relaxed by assuming 
that the velocity be finite at the boundary: The conclusions would be unchanged. We 
want to point out, however, that due to the small viscosity of the electron liquid the 
effective boundary conditions on the electron liquid do not affect considerably its bulk 
motion. 

The analytical stationary solution of Q with these boundary conditions is not 
known. However, one can find an approximate solution by applying the transformation 
x = p, y = -Jkq that maps y = ky/\x\ — 5 into q = \/k(\p\ — 5), q > and transforms 
(jUJ) in a new set of equations where the parameter 1/k appears explicitly. We can 
then expand this equation in powers of 1/k. The zeroth order solution has the form 
v p{Pi o) — 0) v q{Vi q) — ctq 1 ^ 2 {p 2 — S 2 ) where v p ( q ) is the p (q) component of the velocity 
in the p — q plane and a is an arbitrary constant fixed by the requirement that a certain 
amount of charge is flowing through the system. If we transform back to the x — y plane 
we see that the zeroth order solution is given by v x (x, y) = 0, v y (x, y) = a(x 2 — 5 2 ), i.e., 
we have obtained the classical Poiseuille flow [1] . It is well known that the Pouiseille flow 
is stable against small perturbation for almost all Reynolds numbers; The Poiseuille flow 
is laminar up to Re ~ 2 x 10 4 . This flow is linearly stable for any Re, but is unstable for 
non- linear perturbations at large Re: turbulent flow can be observed if, e.g., one could 
force a large enough current. |3J Since for an adiabatic constriction 1/k ~ ( i.e., at any 
given point the system is arbitrarily close to a pipe), we conclude that in an adiabatic 
QPC the flow is laminar for almost any value of Re. We note that compressibility of 
the liquid may instead provide a lower critical Re to observe turbulence. 

Non- adiabatic QPCs - Let us now look at where the constriction is non- 

adiabatic. A simple non-adiabatic potential for which an analytical solution exists is 

% A stability analysis of the solution at large Re numbers shows that the terms of high order appear 
to be important, and one cannot consider only linear perturbations (i.e., proportional to the field 
itself). Assuming non- linear perturbations, one can then show that the Pouiseille flow is unstable for 
an arbitrary (but small) perturbation. This is different from the usual linear perturbation case where 
the system is stable for any small linear perturbation. 
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Figure 1. Critical Re number as a function of the angle j3 of the geometry represented 
in the inset (b). Inset (a) is a schematic of an adiabatic constriction. 



y = k\x\ [see Fig. H](b)]. Microscopically, the point x,y = could be, e.g. a molecule 
sandwiched between two bulk electrodes and current flows from one electrode to the 
other (see e.g. [16J). We are interested here in the dynamics close to (but away from) 
this point. The present problem could also be solved by assuming a finite opening of the 
potential at the origin, i.e. a potential of the form y = k(\x\ —S) (5 > 0). Our conclusions 
would be unaffected by this finite opening. In the case of the system in Fig. [D(b) we know 
there is a critical Re number, R cr , determined by the angle tan(/3/2) = 1/k, above which 
the laminar flow is unstable. A simple calculation gives the implicit relation between 
the angle (3 and R cr [1] 



P(s) = Vl - 2s 2 K(s 2 ) (11) 

= -lz* m + ( 12 , 

o s z s z 

where K and E are elliptic functions and < s < 1 is an arbitrary parameter. This 
critical Re is plotted in Fig. [TJ as a function of the angle (3 and separates the phase 
space in two regions: laminar for Re < R cr , turbulent for Re > R cr . Equation (fl2"j) 
has been again derived with the no-slip condition at the boundary. A finite velocity 
at the boundary (or even a larger velocity at the boundary than at the center of the 
junction) would further reduce R cr for any (3, i.e. it would make the laminar solution 
even more unstable. We can thus conclude that turbulence can be observed in this 
case. Recently the current distribution for a 2DEG QPC has been measured [17]. An 
irregular time-independent pattern has been observed and explained by the presence 
of impurities in the system. We suggest that turbulent effects could be observed in 
similar experiments on the 2DEG if nonadiabaticity is introduced, e.g. by asymmetric 
electrodes that would generate a potential of the type represented in Fig. [0(b). In this 
case we can predict the dimension Ao of the smallest observable eddies in the turbulent 
regime. A simple dimensional analysis gives [U Ao ~ I {R cr /R) m = I {I cr /I) 3/A , where 
I cr is the critical current, i.e., the current which, from (fl2"l) gives R cr , and / the linear size 
of the device. We easily see that Aq decreases rapidly with increasing (3. The viscosity 
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77 can be evaluated through a perturbation theory on the 2D electron liquid [TU]. For a 
two dimensional electron gas (2DEG) in a GaAs heterostructure (m = 0.067m e ) with a 
sheet density of n ~ 10 15 m -2 we have, rj/hn ~ 0.05 [loj"*"! . For a typical current of 1 jj,A, 
Re = 145. We then expect that a turbulent flow is developed for any angle larger than 
(3 ~ it/ 10. For instance, for (3 = rr/2 (R cr ~ 10 from Fig. CD), by assuming a length / of 
a typical device of about 1/iin, we evaluate the typical dimension of the smallest eddies 
to be Ao — 250 nm. Our predictions should be thus readily verifiable experimentally. 

We finally conclude by noting that the geometry of nanoscale structures, like, e.g. 
a molecule between bulk electrodes, is actually closer to a conical structure for which 
the electron liquid is turbulent for relatively small Re [18]. The formation of eddies in 
proximity to an atomic or molecular junction is thus much more likely to occur than 
in the 2DEG case. Moreover, we expect that the finite compressibility of the electron 
liquid will favour the appearance of turbulent flow. 
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